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The Spin-Chern (Cs) was originally introduced on finite samples by imposing spin boundary 
conditions at the edges. This definition lead to confusing and contradictory statements. On one 
hand the original paper by Sheng and collaborators revealed robust properties of Cs against disorder 
and certain deformations of the model and, on the other hand, several people pointed out that Ca 
can change sign under special deformations that keep the bulk Hamiltonian gap open. Because 
of the later findings, the Spin-Chern number was dismissed as a true bulk topological invariant 
and now is viewed as something that describes the edge where the spin boundary conditions are 
imposed. In this paper, we define the Spin-Chern number directly in the thermodynamic limit, 
without using any boundary conditions. We demonstrate its quantization in the presence of strong 
disorder and we argue that Cs is a true bulk topological invariant whose robustness against disorder 
and smooth deformations of the Hamiltonian have important physical consequences. The properties 
of the Spin-Chern number remain valid even when the time reversal invariance is broken. 
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Quantum Spin-Hall insulators represent a new state of 
matter. They were predicted theoretically in Ref. [T] and 
soon after that they were found experimentally.^ Sam- 
ples made of such materials display dissipationless spin 
currents at the edges, that are robust against continuous 
deformations and disorder.^! 

The existence of the edge channels is due to the non- 
trivial topology of the bulk energy bands and two non- 
trivial topological invariants were proposed to describe 
this topology, virtually in the same time: the Z2 invariant 
proposed by Kane and Mel^and the Spin-Chern number 
proposed by Sheng and collaborators^ (first mentioned 
in Ref. In this paper we focus on the later invariant, 
which came under sustained scrutiny because it promised 
a finer classification of the Spin-Hall insulators. This was 
later argued not to be the case.^^ 

The Spin-Chern number was computed by integrat- 
ing the Berry curvature generated by imposing twisted 
boundary conditions on a finite sample. The numerical 
evidence given in Ref. [H] implied that Cs is a robust topo- 
logical invariant. It was later observed, however, that 
one can continuously deform the model using spin ro- 
tations that keep the Hamiltonian's gap unchanged but 
switch the sign of the Spin-Chern number.l^^ This ar- 
gument shows that sometime the structure proposed in 
Ref. E] fails to be a smooth fiber bundle and that Cg may 
not be well defined over the entire Spin-Hall zone of the 
phase diagram. The current understanding is that, when- 
ever one crosses certain zones of the parameter space, Cs 
jumps, but these jumps are always by an even number. 
Therefore, one can still use Cs to formulate a Z2 classifi- 
cation of the Spin-Hall insulators and to efficiently com- 
pute the Z2 invariant. For this reason, the interest in the 
Spin-Chern number continues to be strong. For exam- 
ple, an efficient algorithm for numerical evaluations of Cs 
was proposed by Fukui and Hatsugai,^ and later the al- 
gorithm was used to map Cg for aperiodic systems.^ But 
other works totally dismiss the Spin-Chern number.^^Ei] 



For example Ref. [TT]states that Cs loses its meaning when 
the spin is not conserved. 

In this paper we re-define the Spin-Chern number, this 
time directly for an infinite sample and without involving 
any boundary conditions. This shows, univocally, that 
Cs is a topological invariant describing the bulk and not 
some edge as previous works implied. In our approach, 
the Spin-Chern number is defined as the half difference 
between the Chern numbers for the spin up and spin 
down sectors of the occupied space. We demonstrate 
that the spin up/down sectors can be rigorously defined 
even when Sz is not conserved. 

Furthermore, we state some very general results from 
the non-commuative theory of the Chern numbeil^ and 
then we show how to apply these results to demonstrate 
the quantization of the Spin-Chern number in the pres- 
ence of strong disorder. Several physical consequences 
emerge from this analysis like the existence of critical en- 
ergies below the Fermi level where the localization length 
diverges. Unlike the Z2 invariant, the robust properties 
of the Spin-Chern number remain unchanged when the 
time reversal symmetry is broken. 



I. FIBER BUNDLES AND THEIR 
TOPOLOGICAL CLASSIFICATION 

The classification of the topological insulators is intrin- 
sically related to the classification of the fiber bundles. 
Intuitively, a fiber bundle is a collection of spaces, called 
fibers, that are indexed by a set of parameters that live 
on a manifold. The parameter space is called the base 
manifold and the reunion of all fibers is called the total 
space. An alternative view is to imagine the fibers as be- 
ing attached to the points of the base manifold and then 
glued together by a prescribed topology. 

The fiber bundles encountered in condensed matter 
theory are generated by families of projectors Pe , where 
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0's live on a manifold M like the sphere, torus etc.. 
These projectors act on the same Hilbert space TC and 
they generate the fibers H0=PgH. Since these fibers are 
embedded in the big Hilbert space H, one can tell when 
two fibers He and He' are close to each other and when 
they are not, therefore the bundle of fibers is already 
equipped with a topology. 

Most of the fiber bundles studied in condensed matter 
are generated for finite samples by imposing parametri- 
cally dependent boundary conditions. To be more spe- 
cific, consider a real space 2D lattice with N quantum 
states per site and assume that a tight-binding Hamil- 
tonian H has been defined. The wave functions are N 
dimensional complex vector fields defined over the finite 
lattice: ^(711,712), ni — l,...,Ni and 71,2 — 1,...,-/V2- 
Assume that we diagonalize H by imposing boundary 
conditions like: 

*(7Vi,n2) = Mi(0)*(l,n2) 
*(ni,A^2) = M2(0)*(ni,l), 

where Mi{6) and M2{9) are two NxN matrices depend- 
ing parametrically on 6={9i, 62). Assume that, after the 
diagonalization, we find that the energy spectrum has a 
gap for all the allowed O's. In this case, we can construct 
a fiber bundle by simply considering the projectors Pg 
onto the states with energy below this gap. 

While this is a straightforward construction which per- 
tains to easy numerical implementations, it has a few 
drawbacks which are very often overlooked. The possi- 
bilities of choosing the two matrices seem endless and 
same can be said about the parameter space. How can 
one then discern between a meaningful construction and 
the rest? The ultimate goal of these fiber bundle con- 
structions is to characterize the bulk. Therefore, only a 
rigorous study of the thermodynamic limit can filter out 
the irrelevant matrices and the wrong parameter spaces. 
Most of the constructions will not survive in the ther- 
modynamic limit either because the energy gap closes at 
some point or because the fibers seize to depend smoothly 
on the parameters. Unfortunately, studying the ther- 
modynamic limit requires tedious numerical simulations, 
which most of the time can be carried only up to rela- 
tively small sample sizes. In fact, it is very rare when 
such studies are performed at all. We will discuss in the 
next section how to define fiber bundles directly in the 
thermodynamic limit. 

Wc now turn to the problem of classifying the fiber 
bundles. If the dimension of the fibers Hg is N, the fiber 
bundle generated by a family of projections is a U(7V) 
fiber bundle because two fibers Tig and Tig/ can be linked 
by a unitary transformation Uggi. The question of how 
many topologically non-equivalent U (N) fiber bundles 
can be b uild over a m anifold M has been answered long 
time ago.fi^IIillllinillllt turns out that there is a universal 
manifold X, called classifying space, such that all V{N) 
fiber bundles over M can be generated from continu- 
ous functions / : M X, via the so called pull-back 



construction. Describing this pull-back construction is 
rather technical, but here is a very simple result: two 
fiber bundles over M generated by / and /' can be con- 
tinuously deformed into each other if and only if the two 
functions / and /' can be continuously deformed into 
each other. Therefore, classifying the IJ{N) fiber bun- 
dles over a manifold M is equivalent to classifying the 
continuous functions from M to X . While a fundamen- 
tal result in the fiber bundle theory, this classification 
algorithm is somewhat formal because the space X is 
quite complicate [it is the infinite Grassmannian mani- 
fold GrN{C°°)]. 

A new development in the classification of low dimen- 
sional fiber bundles came from the work of Panati.'^ 
This work provides a computationally more manageable 
classifying approach and states, among other important 
things, that the fiber bundles over 2 dimensional man- 
ifolds without boundaries are completely classified by 
their Chern numbers. Note that the Chern number is 
just one invariant among many topologic invariants that 
can be associated with a fiber bundle. The Chern num- 
ber can be explicitly computed by integrating the Berry 
curvature form: 

dF = Tr{Pg [de, Pg , dg, Pe] A d02 , (2) 

over the base manifold: 

C=^ f dF. (3) 

The result of such integration is always an integer and 
this integer number is called the Chern number. 

Given Panati's result,!^ it is clear that in 2 dimensions 
there are no magic new topological numbers. Therefore, 
the Spin-Chern number must be related to the classical 
Chern number. As we shall see, this is indeed the case, 
but what is new is the structures for which the Chern 
number is computed. 

II. WORKING DIRECTLY IN THE 
THERMODYNAMIC LIMIT 

There is one simple case in which one can work directly 
with infinite samples. This is the case of translational 
invariant Hamiltonians where one can use the Bloch fi- 
bration. The Bloch fibration is a unitary transformation 
from the Hilbert space of the infinite sample into a con- 
tinuum direct sum of N dimensional complex spaces: 

(4) 

where k lives on the Brillouin torus T — [0, 27r] x [0, 27r]. 
U transforms the original Hamiltonian in a direct sum of 
Bloch Hamiltonians: 

UHU-^ = ^H{k). (5) 
feer 
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Obtaining the expHcit expressions of H{k) is, most of the 
time, a very simple exercise. Since U is an isometry, if H 
has an energy gap, then all iJ(A;) have an energy gap. We 
can then compute the projectors P{k) for the states be- 
low this gap and define the fiber bundle ^j^^j- P{k)C^ , 
with the Brillouin torus T as the base manifold. Con- 
sequently, we can compute a Chern number directly for 
the infinite sample. 

Consider now the case when the translational symme- 
try is lost, such as when disorder is present. For sim- 
plicity let us assume that the bulk gap remains open. If 
this is the case, we can define the total space simply as 
PTi, with P being the projector onto the states below 
the gap, but how can we define a base manifold without 
the Brillouin torus? It turns out that the only way out 
of this dilemma is to define a non-commutative Brillouin 
torus, that is a torus whose points are operators instead 
of simple points on a manifold. 

Describing this non-commutative torus is quite techni- 
cal, but we can describe, in quite simple terms, how to 
define the Chern number. Let us show how this defini- 
tion emerges from the translational invariant case. The 
Chern number computed as 

with the Berry curvature given in Eq. |2] (where k takes 
the place of 0) has an expression in the real space: 

C = 27r*tr{P[[ni,P],[n2,P]]}, (7) 

where tr means trace over the quantum states at the 
origin, [, ] denotes the usual commutator and h is the 
position operator: (ni^')(ni, 712) = ni^{ni,n2) {i—l,2)- 
One can see that we don't really need the Brillouin torus 
to define the Chern number! Therefore, one could try to 
define the Chern number for the disordered case using 
this expression. The main question is, will this number 
be an integer? A fundamental result in non-commutative 
geometry says that this is indeed the case after averaging 
over the disorder. 

We would like to state this result in more precise terms. 
For this we need to describe the disorder more precisely. 
To be specific, let us consider the addition of a random 
potential to H: 

K; = ^ A„(w)|n, a)(n, a], (8) 

n.a 

where An(w) is a random variable and a is the index of the 
quantum states (a=l, . . . , N). lo represents a particular 
disordered configuration and can be viewed as a point 
in the disorder configuration space f2, which must be 
equipped with a probability measure dii{ui) to be used for 
averaging over the disorder. We assume that the macro- 
scopic state is translational invariant, in which case the 
lattice translations induce unitary transformations u„ on 
Ti and flows i„ on Q, such that WnK^w^^ = Vt^^j- The 



probability measure dfi{uj) is assumed invariant and er- 
godic with respected to the flows tn . A concrete example 
will be the white noise, for which the disorder configu- 
ration space is the infinite cartesian product of [— |, |] 
intervals: 

x„ez^ (9) 
A point uj of this space is an infinite sequence: 

cj = (. . . ,a;„,cj„+i, . . .) = Xnez2 Wn (10) 

and Aii(cl') is taken as A„(a') = ujn. The probability 
measure d^{ijj) is given by the infinite product of mea- 
sures: dii(uj) — Xngz2 dujti- The ergodic flow is simply 

^m(XneZ2 t^n) — ^neZ'^ ^^n-m- 

With these being said, consider a family of projectors 
P^ acting on the whole Hilbert space TC. We can now 
state the following rigorous result 

Proposition 1. If the matrix element {n\P^\m) decays 
sufficiently fast (exponential decay is enough) with the 
separation \n — m\ and if P^j's satisfy UnPujU^^^ — Pt^uj, 
then: 

C = 2TTi[ dAi(a;) tr{P^[[ni,Pj,[n2,Pj} (11) 

is an integer that is invariant to smooth deformations of 
PJs, as long as they remain localized. 

This is a very general results and it applies to any family 
of projectors and any type of disorder as long as the con- 
ditions inside the Proposition's statement are satisfied. 

If the projectors P^ are associated to the spectral pro- 
jectors onto the occupied states of a Hamiltonian, the 
Chern number defined above gives the disorder averaged 
Hall conductance. Proposition 1 tells us that this av- 
erage remains quantized even in the presence of strong 
disorder. Since the Hall plateaus cannot be explained 
without the strong disorder, it was the above statement 
that completed our understanding of IQHE. 

III. DEFINING THE SPIN-CHERN NUMBER 

The Spin-Chern number was originally defined for a fi- 
nite sample, by choosing Mi{e)^e^^^ and M2{0)^e'^^'^^^' 
in the boundary conditions Eq. [T] Here we define the 
Spin-Chern number directly in the thermodynamic limit. 

Let us start our discussion from the concrete model of 
electrons in graphene;^ 

(mn) ,a 

+iAso E [s • (dfcm X d„fc)]^,^'|m,CT)(n,cr'| 

{{mn)) .aa' 

+iXR J2 [z • (s X dmn)]<7.a' \n, a) (n, a'\. 

{mn) ,(7cr' 

(12) 
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Here, m and n denote the sites of the honeycomb lattice 
and a and a' the electron spin degrees of freedom, tak- 
ing the values ±1. The Hamiltonian acts on the Hilbert 
space H spanned by the orthonormal basis \n,a). The 
simple angular brackets denote the nearest neighbors and 
double angular brackets denote the second nearest neigh- 
bours. In the second sum, k represents the unique com- 
mon nearest-neighbor of m and n. The electrons are 
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considered non-interacting. The three terms in Eq 
are the usual nearest neighbor hopping term, the intrin- 
sic spin-orbit coupling preserving the lattice symmetries 
and the Rashba potential induced by the substrate sup- 
porting the graphene sheet .^^^ We assume that the param- 
eters in the Hamiltonian are chosen so that we are in the 
Spin-Hall zone of the phase diagram.'^ 

There are two sites per unit cell and two spin states per 
site, therefore the model has 4 states per unit cell. Let 
us neglect the disorder for a moment and let P denote 
the spectral projector onto the states below the insulating 
gap of Hq. Under the Bloch fibration U : H ^ ®keT 
this projector becomes: UPU~^=(BkeTP{k). The Bloch 
Hamiltonians H(k) display two upper and two lower 
bands separated by an insulating gap. The Chern num- 
ber of the projectors P{k) onto the lower (and for that 
matter also of the upper) bands is zero, as it will generi- 
cally be for any time reversal invariant band model. Ac- 
cording to Panati's result,'i^the fiber bundle of the occu- 
pied states is trivial, i.e. it has the simplest topological 
structure possible, that of Tx C^. How can there be non- 
trivial topological structures? The answer is because in 
the Spin-Hall effect the non-trivial topological structures 
appear only in the spin sectors, when taken separately. 

We now show how one can define spin up and spin 
down sectors even for the case when Sz does not commute 
with the Hamiltonian. The spin operator is defined as 



|n, cr) = ^<T\n, a) 



(13) 



The goal is to slash the fiber bundle ^^^^PiWjC'^ 
into two non-trvial fiber bundles, that is to make the 
smooth decomposition P(fc)=P_(fc)®P+(fc) for all the 
fc's of the Brillouin torus. At the first sight, this seems 
a very complicated job with no clear chances of suc- 
cess. The key idea is to use the operator Ps^P to do 
just that. Indeed, after the Bloch fibration, UPs^PU^^ 
becomes ®k,zrP{k)szP{k) and we can diagonalize each 
of the operators P{k)szP{k). As we shall see, if the 
Rashba term doesn't exceed a threshold value, the spec- 
trum of P{k)szP{k) consists of two isolated eigenval- 
ues, positioned symmetrically and away from the ori- 
gin, for all the fc's of the Brillouin torus. The spec- 
tral projectors P±{k) onto the negative/positive eigen- 
values are smooth of fc and can be use to achieve the 
decomposition P{k) = P-{k) © P_|_(fc). At this point we 
can can define the Chern numbers C± for the fiber bun- 
dles ®j.g7- P±{k)C^ via Eq.[2]and define the Spin-Chern 
number as C^— ^(C-|_— C_). This is our construction in 
a nutshell. 



We now step back and give a more general construc- 
tion, without involving the Brillouin torus. Assume that 
the disorder is turned on and that we work now with the 
Hamiltonian H(lu)—Ho+XV^. For simplicity, we assume 
that the amplitude of the disorder is not very large so 
that H{uj) still has an energy gap. Let P(w) be the pro- 
jector onto the states with energy below this gap and con- 
sider the operator P{uj)szP{uj)- When the Rashba term 
is zero, Sz commutes with H{oj) and for that reason it also 
commutes with P{uj). Consequently, its spectrum con- 
sists of just two points, ±i. These eigenvalues are highly 
degenerate. When the Rashba interaction is turned on, 
Sz and P(w) no longer commute, therefore the degener- 
acy is lifted and the eigenvalues of P(w)s^P(u;) spread 
towards the origin. As we shall see, if the amplitude of 
the Rashba term does not exceed a threshold value, the 
spectrum of P{lo)szP{lo) remain confined into two iso- 
lated islands, with the origin separating them. Therefore, 
we can define the spectral projectors onto the states with 
eigenvalues below the origin: P_(a;), and above the ori- 
gin: P+{uj). To be able to define the non-commutative 
Chern numbers, we must demonstrate that the matrix 
elements {m\P±{uj)\n) decay exponentially fast with the 
separation |m — n|. 

The exponential decay is required by the conditions 
stated in the Proposition 1. However, we would like to 
pause here and discuss what is the implication of this 
exponential decay for the translational invariant case. 
In this case, we have P±=J P±{k)dk and it is known 
since the work of Kohnis! that only smooth integrals of 
fc lead to exponentially localized functions in real space. 
The reciprocal statement is als o known to be true since 
the work of des Cloizeaux,'2Sllll namely if a total projec- 
tor is exponentially localized then its Bloch decomposi- 
tion is smooth with fc. The smoothness of P±(fc) on k 
is essential because the Berry curvature involves deriva- 
tives over fc. Therefore, the exponential localization of 
{m\P±{ijj)\n) is essential for both translational invariant 
and disordered cases. The good news is that exponential 
localization is always simpler to demonstrate than the 
smoothness with fc. 

Once we demonstrate the exponential localization of 
the kernels, we can define the Chern numbers for the 
spin up and spin down sectors via: 



C±^2t:i / d^i{u:)ir{P±{uJ)[[h^,P±{uJ)],[h2,P±{uJm 
Jn 

(14) 

Since Sz commutes with the translations, we have 



UnP±{uj)u^^ = P±{tr,Uj) 



(15) 



therefore, according to Proposition 1, C± are integers 
that do not change under continuous deformations as 
long as P± (w) remain exponentially localized. As we shall 
see, the exponential localization is guaranteed by the ex- 
istence of the energy gap for H{uj) and by the spectral 
gap of P{uj)szP{uj). To compute the Spin-Chern number 
explicitly, we can turn off the disorder and the Rashba 
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interaction, in which case the model Eq. 12 reduces to 
two decoupled Haldane models!^ for which C±=±l if 
Xso > and C±==Fl if A50 < 0. 

We now start the proof of the exponential localization 
of the projectors P±{uj). To ease the notation we drop 
uj. First, let us show that the gap Ai of the operator 
PszP, viewed as an operator on the space of occupied 
states K.=PTL, remains open and clean when the Rashba 
term is turned on, as opposed to immediately closing or 
filling with additional spectrum due to some instability. 
For this, we notice first that the projector P in the occu- 
pied space is analytic in A/j. This property is protected 
by the insulating gap A2 of the Hamiltonian H{uj). We 
then show that, at least for small Xn, Ps^P ~ <^Ik is in- 
vertible for C in the vicinity of [here Iiq is the identity 
operator for the space /C]. One could try to work with 
the expression 



{Ps^p - qik.)-^ - p{Ps,p - ar 



(16) 



and use the continuity of P, but this expression has a 
problem when C,—0 since {PszP — CI)~^ diverges because 
we include the un-occupied states where PszP is zero. 

Here is an alternative approach inspired from Ref. 1231 
Let i?(C) = P{sz — C)~^P, with ^ in a neighborhood of 
zero. We have successively: 



P{Sz 



R{C){PszP-CIk)^ 

-0-'{{Sz~C)P+[P,sz]}P 



= p + P{h - 0-^[P, h]P ^Ik + Q. 

Q is small, at least for small A/j, since 

Q = P{sz^O-\P-Po.Sz]P. 



(17) 



(18) 



where Po is the projector onto the occupied states for 
Ai? = 0. The small factor comes from P — Pq, which 
is proportional to A^. Above, we used the fact that 
commutes with Pq, even in the presence of disorder. In 
this case, the operator Ijc + Q has an inverse given by 
the convergent series J2^=o(~Q)" ^^"-^ find: 

(/^ + Q)-'R{0{PszP - CIk) = Ik, or (19) 



C- 
— <- 



spectrum of PszP 



-1/2 



Ai/2 



spectrum of PszP 
^ 



1/2 



Ret; 



FIG. 1; Illustration of the spectrum of Pa^P for small Rashba 
coupling, of the gap Ai and of the integration contours C± 
used in the main text. 



then the following is true.lSl! 

Proposition 2. If T is exponentially localized, then 



UqTU-q is a bounded operator for all orientations 



of q, provided |q| is smaller than a, and conversely: if 
Tq is bounded for any orientation of q and \q\ < a, then 
T is exponentially localized with a localization exponent 
equal or larger than a. 

Also, it is a fact that if T is exponentially localized, then 
[II II denotes the operator norm] 



lir-r„ 



as q ^ 0, 



(23) 



i.e. the difference between T and Tq is small for q small. 
Since 



P±^^ f P{PszP-Ci)-'dC 

Jc± 



(24) 



it is enough to show the exponential localization of 
P{PszP — CI)~^ for all C on the contours C± shown in 
Fig. [1] Now the projector P itself is exponentially local- 
ized and a general proof of this statement can be found 
in Ref. 25. Then Pq is bounded and we have successively: 



UqP{PSzP - Ciy'Uq = PqiPqSzPq " C^) " 



and 



PgSzPq -CI = PSzP - CI 



(25) 



(26) 



{PszP - (Ik)-' = {lie + Q)-'P{sz - Cr^P (20) 

Thus we showed that resolvent of PszP is finite for ^ near 
the origin, which excludes the existence of any eigenval- 
ues in this region. 

We now show that, as long as the gaps Ai and A 2 of 
PszP and II(lu) remain opened, the spectral projectors 
P± are exponentially localized. By exponential localiza- 
tion of an operator T we mean the existence of a strictly 
positive a such that 



\{m,a\T\n,a') \ < ct. e""!'"-"!. 



(21) 



thus 



P H P 



ail > ^ - 

S II - 2 



PqSzPq - PSzP\\ 



(27) 



which is strictly positive for small q due to Eq. 23 In 
other words, PqS^Pq — C,I is invertible, at least for small 
q, which shows that all operators appearing in the right 



hand side of Eq. 25 are bounded. Consequently, P± {u) 



are exponentially localized. 



IV. DISCUSSION 



We will use the following simple observation: if Uq de- 
notes the non-unitary transformation 



Uq\n,a) = e'?'"|n,CT) 



(22) 



Our results show that the Spin-Chern number is a 
topological invariant that is well defined in the thermo- 
dynamic limit. It is now clear that the robustness is 
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due to the existence of two spectral gaps: the insulat- 
ing gap of the Hamiltonian and the spectral gap of the 
operator Ps^P. Contrary to many folkloristic believes, 
not all topological invariants are robust against disorder. 
In fact, to date, a rigorous proof for such statement ex- 
ists only for the classic Chern number. The proof makes 
heavy use of the particular expression of the Chern num- 
ber, thus it cannot be automatically extended to other 
topological invariants. Therefore, it is remarkable that 
we can add the Spin-Chern number to the list of topo- 
logical invariants that are robust to disorder. 

The robustness of the Spin-Chern number has physi- 
cal consequences. We recall that, according to Ref. HHI 
and I27|, generically for dimension smaller than two, all 
states are exponentially localized at any strength of dis- 
order; for dimension greater than two, there exist ex- 
tended states for low disorder and for dimension two all 
the states are localized except for states corresponding to 
isolated critical energies at which the localization length 
diverges. The Chern number in IQHE not only has a 
classifying role, but its quantization in the presence of 
disorder implies the existence of such critical energy re- 
gions below the Fermi level.l^ The quantization of the 
Spin-Chern number leads to the same conclusion, which 
explains the finite Spin-Hall conductance in the presence 
of strong disorder."^ This remarkable property is solely 
related to the spin, because the Spin-Chern number is 
robust against perturbations that violate the the time 
reversal symmetry. Indeed, none of our constructions or 
proofs use the inversion symmetry as an argument. 

The operator is not special. In fact, Cs remain in- 
variant to deformations of the Hamiltonian and of Ps^P, 
as long as the spectral gaps of the two operators remain 
open. Closing any of the two gaps can result in jumps 
for the the Spin-Chern number. The jumps are always by 
an even number. The argument against the Spin-Chern 
number in Refs. |7] and H] was that, using rotations in 
the spin sector, one can connect the Hamiltonians with 
+\so E^nd with —XsOi without changing the insulating 
gap. But at the end of such rotation, Cg changes sign 
and it was concluded that Cg is not a well defined topo- 
logical invariant. The change of sign was attributed to 
the closing of the gap of the Hamiltonian with the spin 
boundary conditions. We now can give an alternative 
explanation: during the rotation in the spin sector, the 
gap of PszP closes and then opens again. This problem 
can be easily fixed, the solution being to deform not only 
the Hamiltonian but also the operator PszP- For the 
case of spin rotations, this can be easily accomplished by 
applying the spin rotations to Sz inside PszP- If there is 
such easy fix, then the old question resurfaces: does the 
Spin-Chern number contain more information than the 
Z2 invariant? The answer is no. After the continuous ro- 
tation mentioned above, PszP — > —PszP and P± Pzp. 
Thus, there is no canonical way to chose the projectorss 
P± and, since this choice determines the sign of Cg, the 
sign contains no additional information. 

The construction can be easily generalized. For more 



complex insulators, it is very probable that one can re- 
place Sz by other non-trivial operators w by combining 
spin and/or point group symmetry operators. Using such 
operators one might discover nontrivial topological struc- 
tures in seemingly trivial insulators. For example, it was 
recently pointed out^- that certain surface states in or- 
dinary semiconductors can have topological origins. The 
non-zero Chern numbers for the spectral sectors of PwP 
will then lead to a t&-Hall effect. For other models, we 
can find operators PwP that have n islands of isolated 
spectrum with exponentially localized projectors. If w 
commutes with the translations of the unit cell, we can 
define a Chern number for each spectral island of PwP. 
A simple example of this type is the model of Eq. [12] with 
spin I particles. For Afl=0, Sz commutes with the Hamil- 
tonian and PTi splits into four sectors, corresponding to 
Sz — —\, \ and I . Turning on the disorder will not 
affect these eigenvalues. For each sector we can define a 
Cern number, which take the values (assuming A50 > 0)^ 

C_3 = C_i = -1 and Ci = C3 = 1. When the Rashba 

22 22 

interaction is turned on, the spectrum of PszP spreads, 
but is still contained into four isolated spectral islands. 
Thus, we can still split PTi into four sectors and our anal- 
ysis shows that the Chern numbers for each sectors will 
be conserved as long as the the insulating gap remains 
open and the spectral islands of PszP remain isolated. 
The spin-Chern number can be defined in various ways, 
depending how we group the sectors. If we repeat the 
construction for spin ^ and put the negative Sz sectors 
into /C_ and the positive Sz sectors into /C^., Cs becomes 
2. Other possible groupings gives Cs = and Cg — —2. 

The existence of the nontrivial Chern numbers for dif- 
ferent sectors does not automatically imply the existence 
of chiral edge modes. Generically, only the insulators 
with odd Cs display such edge modes. The discussion 
given here hardly touches the problem of classification 
and of the edge modes. For these two problems, one has 
to explore how the occupied states relate (connect) to 
the un-occupied states. This has to be done on solvable 
models. It is at this point where our analysis becomes 
relevant because now we have a guiding principle which 
tells when the realistic models, which should include dis- 
order and will generally not be solvable, can be deformed 
into smooth solvable models without changing the topo- 
logical invariants. The robustness of the edge modes was 
also recently investigated in Ref. '29! which introduced a 
quantized edge index. Although this edge index was also 
constructed via a splitting, the connection between the 
edge index and Cs is not clear to us at this moment. 

We want to make a final remark about the edge modes. 
The existence of nontrivial topological sectors can indi- 
cate reach edge and surface physics. For example, for the 
model of spin | particles, when cutting an edge we can 
be sure that four [which become 2] edge bands shoot out 
of the bulk spectrum. The bands hybridize and return 
back into the same part of the bulk spectrum were they 
originated. Nevertheless, these bands lead to [gapped] 
edge states, which can still be useful for practical appli- 
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cations. In the presence of disorder, these edge bands will 
localize, but if the localization length is large, these edge 
states, for example, can efficiently trap light and thus be 
useful in photovoltaic devices. 

V. CONCLUSIONS 

We have demonstrated that the Spin-Chern number 
can be defined in the thermodynamic limit without the 
need of any boundary conditions, therefore showing that 
the Spin-Chern number describes the bulk and not the 
edges as it is currently believed. By making a connec- 



tion with the non-commutative theory of Chern number, 
we were able to demonstrate the robustness of the Spin- 
Chern number against disorder and smooth deformations 
of the models. The robustness, which is solely related to 
the spin and has nothing to do with the time reversal 
symmetry, has physical consequences: it implies the ex- 
istence of critical energies below the Fermi level where 
the localization length diverges. We have also theorized 
possible generalizations of our construction, which could 
aid the search and discovery of new topological phases. 
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